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Abstract 

We analyze the mean-squared error (MSE) performance of widely linear (WL) and conventional subspace-based 
channel estimation for single-input multiple-output (SIMO) flat-fading channels employing binary phase-shift-keying 
(BPSK) modulation when the covariance matrix is estimated using a finite number of samples. The conventional 
estimator suffers from a phase ambiguity that reduces to a sign ambiguity for the WL estimator. We derive closed- 
form expressions for the MSE of the two estimators under four different ambiguity resolution scenarios. The first 
scenario is optimal resolution, which minimizes the Euclidean distance between the channel estimate and the actual 
channel. The second scenario assumes that a randomly chosen coefficient of the actual channel is known and the 
third assumes that the one with the largest magnitude is known. The fourth scenario is the more realistic case 
where pilot symbols are used to resolve the ambiguities. Our work demonstrates that there is a strong relationship 
between the accuracy of ambiguity resolution and the relative performance of WL and conventional subspace-based 

, , estimators, and shows that the less information available about the actual channel for ambiguity resolution, or the 

f-H lower the accuracy of this information, the higher the performance gap in favor of the WL estimator. 

Index Terms -Widely Linear, Subspace, SIMO, Channel Estimation. 



I. Introduction 



Subspace-based estimation is one of the most popular approaches for blind channel estimation. Originally 
■^j- proposed in [ 1 ], subspace-based methods estimate the unknown channel by exploiting the orthogonality between the 
^ signal and noise subspaces in the covariance matrix of the received signal, offering a convenient tradeoff between 
performance and computational complexity (2j. They potentially outperform methods that are based on higher-order 
^ statistics when the number of available received signal samples is limited, since second-order statistics (SOS) can 
q be estimated more robustly in such conditions [3J. Subspace-based channel estimation has been applied in a wide 
variety of communication systems, including single-carrier systems JTJ, multicarrier (OFDM) systems Q, (5J, and 
T 7"! spread spectrum (CDMA) systems (6j-|[8j. 

It has become a well-known fact that widely linear (WL) processing (9J, which operates on both the received 
^ signal and its complex conjugate, can improve the performance of SOS-based algorithms when the signal under 
consideration is improper. Improper signals, which are characterized by a non-zero pseudo-covariance, result when 
communication systems employ real modulation schemes such as amplitude-shift-keying (ASK), binary phase-shift- 
keying (BPSK), minimum-shift-keying (MSK) and Gaussian minimum- shift-keying (GMSK). By augmenting the 
observation space, WL processing is able to access and exploit the information in the pseudo-covariance of the 
signal. This has prompted researchers to propose WL versions of subspace-based channel estimation algorithms in 
order to improve channel estimation accuracy in communication systems that employ real signaling. In the context 



of DS-CDMA systems, a WL subspace-based channel estimation algorithm was proposed in 1 1 1 for the case of 
BPSK modulation. In addition to exhibiting superior mean-squared error (MSE) performance, the WL algorithm 
was able accommodate almost twice as many users. Similar observations were made in the context of multicarrier 
CDMA in fill . A WL subspace-based algorithm was developed for the estimation of single-input single-output 
(SISO) FIR channels with improper input signals in (12}, showing superior mean-squared error (MSE) performance 
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to the conventional subspace-based method implemented using an oversampling factor of 2. WL subspace-based 
channel estimation was also used in the context of interference-contaminated OFDM systems in |l3j. 

Conventional subspace-based channel estimation suffers from an inherent phase ambiguity which can only be 
resolved through the use of additional side information. An important advantage of WL estimation is that it reduces 
the phase ambiguity into a sign ambiguity |13|, |14|, which is intuitively easier to resolve. For a meaningful 
comparison of WL and conventional channel estimation, both the phase ambiguity and the sign ambiguity have 
to be resolved. Unfortunately, there is no agreement in the literature on how to resolve the two ambiguities. In 
many works on blind channel estimation, it is assumed that the one of the channel coefficients, typically the first 
one, is known Q, (15), p6*[ . This assumption was used to resolve the ambiguities of the conventional and WL 
estimator in fOJ. In (12), the channel coefficient with the largest magnitude was assumed known. Other works on 
WL subspace-based estimation fTOj , fl4j do not explicitly mention how they resolve the phase and sign ambiguities. 

In this work, we consider the problem of blind estimation of single-input multiple-output (SIMO) flat-fading 
channels [ 17 1, [ 18 ]. We assume that BPSK modulation is used, resulting in an improper received signal. Our goal is 
to examine whether and under what conditions WL subspace-channel estimation constitutes an appealing alternative 
to conventional approaches. We pay special attention to the practical situation where estimation is performed with a 
finite number of received samples, and we derive highly accurate closed-form expressions for the MSE performance 
of the conventional and WL channel estimation algorithms. Our MSE analysis explicitly takes into account the effects 
of phase and sign ambiguity resolution by considering four different approaches to resolving the said ambiguities. 
The first approach we consider is optimal phase and sign ambiguity resolution where the applied phase and sign 
correction minimizes the Euclidean distance between the channel estimate and the actual channel. For this case, we 
also derive a closed-form expression for the probability that the WL estimator outperforms the conventional one, 
when the statistics of the channel are taken into account. The second approach assumes that a randomly chosen 
channel coefficient is known, and the third approach assumes that the channel coefficient with the largest magnitude 
is known. For the latter, we also derive a lower bound on the probability that the WL estimator outperforms the 
conventional one. Interestingly, the relative performance of the two estimators is different depending on which of 
the above three approaches is adopted. Under optimal correction, the conventional estimator slightly outperforms 
the WL estimator. However, the WL estimator is significantly better than the conventional estimator when the first 
channel coefficient is assumed known, and slightly better when the channel coefficient with the largest magnitude is 
assumed known. All these approaches assume that certain information about the actual channel is perfectly known 
at the receiver. In practice, however, pilot symbols have to be used to estimate such information, and the incurred 
estimation error will contribute further to the resulting MSE. We therefore consider a fourth approach in which 
a limited number of pilot symbols is used to resolve the ambiguities of the two estimators, and we derive the 
corresponding MSE expressions. As it turns out, the WL estimator is significantly better, coming very close to 
optimal performance even when a single pilot symbol is used. This conforms with the intuition that sign ambiguity 
is much easier to correct than phase ambiguity. The four scenarios demonstrate that the less information available 
about the actual channel for ambiguity resolution, or the less accurate this information is, the higher the performance 
gap becomes in favor of the WL estimator. 

Our work is the first to present a thorough analytical study on the relative MSE performance of conventional 
and WL subspace-based estimators. To the best of our knowledge, we are the first to observe and analyze the 
close relationship between the accuracy of phase and sign ambiguity resolution and the relative performance of 
WL and conventional subspace-based estimators and to analyze theoretically the practical scenario when pilots are 
used to resolve the phase and sign ambiguities. As such, our work offers new and unique insights into the relative 
performance of conventional and WL subspace-based channel estimation algorithms. 

The remainder of the paper is organized as follows. In Section [TTJ we introduce our system model and the 
conventional subspace-based channel estimation algorithm for the SIMO flat-fading channel model. The WL 
subspace-based channel estimation algorithm is developed in Section III The MSE performance of the conventional 
estimator under finite sample size is derived in Section IV for the four scenarios of phase ambiguity resolution. We 
derive the MSE performance of the WL estimator for the corresponding four sign ambiguity resolution scenarios 
in Section [Vj In Section VI we derive the probability that the WL estimator outperforms the conventional under 
optimal ambiguity resolution, as well as a lower bound on the probability that the WL estimator outperforms 
the conventional one when the channel coefficient with the largest magnitude is known. In Section VII| we use 
Monte-Carlo simulations to verify the accuracy of our analytical results and to compare the performance of the 
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Fig. 1. The single user SIMO Communications Model. 



two estimators in the four different cases. Finally, we present our conclusions in Section VIII 



II. System Model and Background 



We consider the single user SIMO flat-fading model used in [ 18 1-|20] and illustrated in Fig. [T] In the zfh symbol 



period, the transmitter sends a binary -phase-shift-keying (BPSK) data symbol b(i) taking the values ±1 with equal 
probability. For a J-antenna receiver, the corresponding J x 1 received vector is 

r(i) = b{£)g + n{i) = b{i)\\g\\h + n(»), (1) 

where g = [g\, . . . ,gj] T is the unknown vector of fading channel coefficients, h = g/\\g\\ = [hi, . . . , hj] is the 
normalized (to unit-norm) version of the channel, and n(i) is the complex additive white Gaussian noise (AWGN) 
vector with mean zero and covariance K[n(i)n(i) H ] = a 2 1. The fading coefficients gi, . . . , gj are modelled as 
independent and identically distributed (i.i.d.) CM(0, 7 2 |j]and are assumed to remain fixed throughout the estimation 
period. The transmit SNR in dB is —10 log(o" 2 ). It is well known that SOS-based algorithms are blind to the channel 
magnitude \\g\\, i.e., they can only estimate the direction and not the norm of g. For this reason, we will focus 
on estimating h, and, with a slight abuse of terminology, we will also refer to the elements of vector h as fading 
coefficients. 

Let R = K{r(i)r(i) H } be the covariance matrix of the received signal, then 

R = gg H +a 2 I = \\g\\ 2 hh H + a 2 1 . (2) 

As we can see from ((2]), the matrix R has two distinct eigenvalues, Ai = (||gf|| 2 + <r 2 ) and A2 = cr 2 , the latter having 
a multiplicity of J — 1. The normalized channel h is an eigenvector of R corresponding to the largest eigenvalue. 
If we let it be a unit-norm eigenvector of R corresponding to Ai, then u has the form u = he^, where j = y/—l 
and <p £ [0, 2ir) is an arbitrary angle which represents the ambiguity in the phase of the channel estimate which 
is common to all SOS-based estimators. In practice, R is not known beforehand and is commonly replaced by 
its sample-average estimate R = 4 Y2i=i r(i)r(i) H , using a finite sample of ./V received vectors, r(l), . . . , r(N). 
In this case, the channel is estimated by the unit-norm eigenvector corresponding to the largest eigenvalue of R. 
We denote by u = [ui, . . . uj] T the subspace-based channel estimate prior to phase correction. For a meaningful 
study of the performance of the subspace-based estimator, it is essential to resolve the phase ambiguity by applying 
an appropriate phase shift to u. The most common convention is to assume that one of the channel coefficients, 
typically the first one, is known (5J, JT3J, fT6| ], and rotate the vector u such that the corresponding component has 
the same phase. This approach may lead to inaccurate results when the chosen coefficient happens to have a very 



small magnitude [21 1. An alternative approach which leads to more accurate ambiguity resolution is to assume that 
the known channel coefficient is the one with the largest magnitude (12). In both cases, however, the phase shift is 
suboptimal because it uses information from only one channel component. In fact, the optimal phase shift which 
minimizes the Euclidean distance between the channel estimate and the actual channel is the phase of their inner 



'The notation CAf(/i, r ) 2 ) is used to refer to the complex Gaussian distribution with mean /i and variance r y 2 . 
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product pT|. In practice, however, none of the information used in these three approaches is readily available at 



the receiver, and pilot symbols have to be used to resolve the phase ambiguity. Obviously, the choice of how to 
resolve the phase ambiguity will affect the observed mean-squared error (MSE) performance. When training pilots 
are used, the error in estimating the desired phase shift will further contribute to the MSE. 

In section |IV| we will derive accurate closed-from expressions for the MSE error of the conventional subspace- 
based estimator under each of the four scenarios. 

III. The Widely Linear Subspace-based Estimator 

Due to the use of real (BPSK) modulation, the received signal is improper, and widely linear processing may be 
applied to access the information in the pseudo-covariance of the signal. In this section, we will develop the widely 
linear version of the subspace-based channel estimator described in Section [n] Widely linear processing operates 
on an augmented vector r(i) formed by stacking the received vector r(i) and its complex conjugate r(i)*. The 
vector r{i) is given by 



r(i) 



r{i 
r(i) 



b(i)\\9\ 



where h = [ h 2 



h H ] T and h(i) 



h 

h* 

[n(i) T n{i) H ] T 



+ 



n(i) 



nil 



b(i)\\g\\h + n(i) t 



R 4 E[r(i)r(i) H 



The co variance matrix R of r(i) is given by 

R C 

C* R* 



(3) 



(4) 



where C = K{r(i)r(i) T } is the pseudo-covariance of r(i). It is easy to see that the augmented channel vector h 
is an eigenvector of R corresponding to the largest eigenvalue, X u , = 2\\g\\ 2 + a 2 . 

WL channel estimation reduces the inherent phase ambiguity into a sign ambiguity. This fact becomes clear 
by reformulating the WL subspace-based estimator into the following equivalent real representation. Let f(i) = 
[$l{r(i)} T , Q{r(i)} T ] T be the real representation of the received vector r(i), and h = [$l{h} T ', ^s{h} T ] T = 
[hi, ... , Ji2j] T be the real representation of the channel h. Finally, let 



Then, h = tyh and r{i) = S&f(i), or, equivalently, h 

f(i) is R = E{r(i)f(i) T } 



Ij 
Ij 



Jlj 
-Jlj 



\^ H h and 



00 



2 * 



g\\ hh + \I, and it is related to R by R 



(5) 

rU). The covariance matrix of 
Moreover, there is a 



1-1 correspondence between the eigenvalues and eigenvectors of R and those of R. Hence, WL subspace-based 
channel estimation can be performed in the real domain through the eigendecomposition of R. We are left with a 
sign ambiguity because both h and — h are unit-norm eigenvectors of R corresponding to its largest eigenvalue. 
For the remainder of the paper, we will use the real representation of the WL subspace-based channel estimator. 

Since R is not known to the receiver beforehand, it is approximated by its sample-average estimate R = 
Jf Y^iLi r(i)r(i) T . To estimate the channel, we obtain the unit-norm eigenvector ti corresponding to the largest 
eigenvalue of R. This channel estimate suffers from a sign ambiguity which requires extra information to be 
resolved. In |T2j , the component of largest magnitude in the real representation of the channel was assumed to 
be known and used to resolve the ambiguity, while in fOJ, the real representation of the channel was converted 
back into a complex vector and the phase of the first coefficient of the complex channel was assumed known and a 
phase shift was applied to match the phase of the two vectors. The approach we follow in our work is to perform 
sign-ambiguity correction directly to the real channel estimate. We maintain fairness in the correction of the two 
ambiguities (phase in the case of conventional and sign in the case of WL) by assuming that the same information 
is available to both estimators. Hence, we will we perform sign ambiguity resolution under the same four scenarios 
discussed in Section [IT] Closed-form expressions for the MSE in all four cases for a finite number of samples will 
be derived in Section |V] 

IV. MSE Performance of the Conventional Estimator 

In this section, we derive closed-form expressions for the mean-squared error performance of the conventional 
subspace-based estimation under the four assumptions on phase ambiguity resolution that were briefly discussed 
in Section [n] It is more convenient to start our analysis with the case of optimal phase correction, as the resulting 
MSE expression will be used in the derivation of MSE expressions for the other three cases. 
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A. Optimal Phase Correction 

Let 9 D be the optimal phase shift that minimizes the Euclidean distance between the channel estimate and the 
true channel. It is straightforward to check that 9 Q is given by 

6» = argmin lle^u - h\\ 2 = Z (u H h) . (6) 

6»e(0,27r] V ' 

Let h Q = e 3 "ii = [h 0t i, . . . , h Q j] be the optimally phase-corrected estimate of the channel, and let 5h = h Q — h 
be the corresponding error in the estimation of h under optimal phase correction. We can decompose 5h as 

5h Q = q Q + a h, (7) 

where q Q = [q ,i, ■ ■ ■ , qo,j] T is the component of 5h a that is orthogonal to h and a Q is a complex scalar. Hence, 

\\5h \\ 2 = \\q \\ 2 + \a Q \ 2 . (8) 

A closed-form expression for E {||S/i 2 ||} is given by the following theorem. The proof can be found in Appendix A 



Theorem 1 Under optimal phase correction, \cxq\ 2 <C ||<? || 2 and 

E{||^ || 2 }^E{||g || 2 }^^l^(a 2 || ff || 2 + f 7 4 )(J-l). 

As we shall see shortly, the MSE under optimal phase correction lower bounds the MSE for the other three scenarios. 
Since the term a Q h has a negligible impact on the estimation error 5h a , we will assume for the remainder of our 
work that a Q = 0. 



B. Suboptimal Phase Correction 

The phase ambiguity can be resolved suboptimally if only one of the channel coefficients is known. Let he, £ £ 
{1, . . . , J}, be the known channel coefficient. The suboptimally-corrected channel estimate is h s = ue 303 , where 
9 S = Z (u}hi). The two channel estimates h s and h Q are related by h s = h e^ 6s ' 6o \ and the estimation error is 
given by 

5h s = h s -h = h o e? ( - '- ') - h 

The resulting MSE is now 

E{||^ s || 2 }=E{|| Qo || 2 }+E{|e^)-l| 2 } 

= E {||g || 2 } + 2 - 2E {cos (9 s - 9 )} . 

The closed- form expression for E{||q || 2 } is available in Theorem [TJ so it remains to find a closed-form for 

E {cos (9 S — 9 Q )}. Since v,£ = [hg + q ^)e~ 3d °, we obtain 

9 S - 9 D = Z (u* e h e ) - Z (u H h) = Z (q* o/ h e + \h e \ 2 ) , (11) 
It is shown in Appendix [A] that 

q ~ -l—VV H 6Rh, (12) 
llffll 

where V is a J x (J - 1) matrix whose columns are the orthonormal eigenvectors of R corresponding to the 
eigenvalue a 2 , and that 

MM 2 } * j^f^ 2 \\9\\ 2 + ^ 4 )(i - M 2 )- (13) 



Using (12 1, it can be shown that ~E{q 0i e} ~ and that q Qy £ is the sum of N independent and identically distributed 
(i.i.d.) random variables. Furthermore, it is also shown in Appendix B 



that E jg 2 ^j ~ 0. Invoking the Central 



Limit Theorem, we can approximate q £ by a (proper) complex Gaussian random variable with the same mean 
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and variance. Under this approximation, the term q*pht + |/i^| 2 is also a proper complex Gaussian random variable 
with mean \h(\ 2 and variance 1 4 (cT 2 ||ff|| 2 + a 4 )(\he\ 2 — |/if| 4 ), and furthermore, i? = Z ylg^hi + |/if| 2 ^ has 
the Ricean phase distribution [22| with the following probability density function 

= -L e -" (l + ^pcos(7?)e pcos2W [l + erf(pcos(0))]) , 



2vr 



where 



(a2|| ff p + a 4 )(1 _|^| 2) 



and erf(-) is the error function [23]. Then, E{cos(?9)} is given by [24| 



E{cos(i?)} 




e 2 



/n 



(14) 
(15) 

(16) 



where Iq(-) and Ii(-) are Modified Bessel Functions of the First Kind of orders zero and one, respectively |23|. 
Going back to (TTOb, we obtain the closed-form MSE expression 



E{\\5h s \\ 2 } 



N\\g 4 



-(a 2 \\g\\ 2 + a 4 )(J-l) + 2 




e ^ 



in 



+ h 



This expression can be simplified if we use the following approximation [25] 

E{cos(tf)} ~ e - ^ = e's^F^^H'+^KiSii* -1 ), 

that yields 

E { \\5h s \\ 2 } ~ 1 ^-^{a 2 \\g\\ 2 + a 4 )(J - 1) + 2 - 2e -^ (<T2||fl||2+<T4) (^- 1 ) 
C. Largest-magnitude Phase Correction 



(17) 

(18) 
(19) 



We can see from eq. ( [19] ) that the MSE decreases as the magnitude of the fading coefficient hi increases. For a 
fixed channel g, the MSE would be minimized when the fading coefficient with the largest magnitude is known. We 
let L be the index of the channel coefficient with the largest magnitude, and denote by h a and 5h a = h a — h the 
resulting channel estimate when hi, is known, and the corresponding estimation error, respectively. The resulting 
MSE is thus 

1 



MH^all 2 } 



2 \\g\\ 2 + a 4 )(J-1) + 2- 2e 4N " 



(20) 



We can obtain a simpler closed-form for E {||<5h. a || 2 } by using the Taylor series expansion of the exponential term 
in (20 1. Since |/il| 2 > !/«/» we have that 



4JV||flf| 



(^llfll| 2 + ^ 4 ) 



hr\ 2 



1 < 



J - I 



a 



N\\g\\ 2 N\\g\ 



(21) 



The RHS of (21) is typically much smaller than 1, so we can accurately approximate E { ||5ft a || 2 } by using the 
first two terms in the Taylor series expansion of the exponential term, resulting in the closed-form 



E{\\5h a \\ 2 } ~ 



N\\g\ 



: (o- 2 \\g\\ 2 + a*) J + 



2\h L \ 2 2 



(22) 



We will use eq. ([22]) later on to obtain bounds on the probability that the WL estimator outperforms the conventional 
estimator under largest-magnitude phase and sign correction. 
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D. Training-based Phase Correction 

The previous three conventions for phase correction assume that certain information about the channel is perfectly 
known. In practice, however, no such information is available and pilot symbols have to be used to resolve the 
phase ambiguity. We will now investigate the resulting MSE if K pilot symbols are used to estimate the optimal 
phase shift O , where the transmitted symbols are set to +1. Let z^, k = 1, . . . , K , be the received sample vectors 
during the training period. Thus 

Zk = 9 + n k , (23) 

where rik,k = 1,...,K, are the complex AWGN vectors with mean zero and covariance a 2 1. To reduce the 
variance of the noise, we average the received vectors during the training period, obtaining 

K 

Z m = 9 + T7 > n k- (24) 



k=l 



Using z m , we estimate 9 Q by O = Z (u H z m ), and denote by e = 6 



the error in estimating O . We denote 



by h t = Me j9 ° = h e J£ and 5h t = h t — h = q e JE + (e J£ — l)h the phase-corrected channel estimate, and the 
corresponding estimation error respectively. The resulting MSE is given by 



E{||(5^|| 2 } =E{||q || 2 } + 2-2E{cos(e)}. 
Moreover, since u = e~ 3do h , we obtain 

H 



ii H z 



°j e °h: 



9 + 



~ J2 n k) = e36 ° [ HSU + J( kH o n *) ' 
k=l / \ k=l / 



(25) 



(26) 



which means that e = ^(11^11 + -kzZk=i^o n k)- The quantity \\g\\ + zZk=i n k is a (proper) complex 



Gaussian random variable with a mean of \\g\\ and a variance of a /K, and the error e has the Ricean phase 



distribution |22| with parameter /3 = 

E{cos(e)} 



^llsll 



. Hence, 



7T/3 




+ h 



r^J g 4/3 



(27) 



and 



E{||^|| 2 } 



N\\g\ 



(a 2 \\g\\ 2 + a 4 )(J-l) + 2 



irK\\g\ 



KM 
2a 2 



+ h 



K\\g\\ 
2a 2 



(28) 



1 



l 2 1 1 ||2 



+ cx 4 )(J-l) + 2-2e 



N\\g\\ A 

As expected, the MSE under training-based phase correction approaches the one under optimal correction as the 
number of pilots K increases. 



V. MSE Performance of the Widely Linear Estimator 

In this section, we will analyze the MSE performance of the WL subspace -based estimator, where the phase 
ambiguity is replaced by a sign ambiguity. In resolving the sign ambiguity, we will consider the same four scenarios 
that were considered in resolving the phase ambiguity of the conventional subspace-based estimator. Providing the 
same information to the WL estimator as was made available to the conventional estimator guarantees a fair 
comparison of the performances of the two estimators. 
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A. Optimal Sign Correction 

Let b Q be the optimal sign correction that minimizes the Euclidean distance between the channel estimate and 
the true channel. It is straightforward to check that 

b Q = arg min libit — h\\ 2 = sgn (ti T h) . (29) 

be{±i} x ' 

We also note that b Q = sgn {cos 6 }, where O is the phase shift we used to do optimal phase correction. Denoting 
by h = b Q u = [h 0t i, . . . h ,2j] an d 6h a = h a — h the optimally sign-corrected channel estimate and the resulting 
estimation error under optimal sign correction, respectively, the error 8h can be decomposed as 

5h a = q Q + p a h, (30) 

where q Q = [q ,i, ■ ■ ■ , q ,2j] is orthogonal to h and p Q is a real scalar. Hence, ||5h. || 2 = ||q || 2 + | /x 1 2 . The 
following theorem presents a closed-form expression for the mean-squared error E |||Jh, || 2 }. 

Theorem 2 Under optimal sign correction, |/i G | 2 -C ||q I| 2 > an d 

V{\\5h f} ^E{\\ Qo f} * (yll^l| 2 + (2J-1). 

The proof of the above theorem is found in Appendixjc] Since the term fj, a h has a negligible impact on the estimation 
error 5h , we will assume for the remainder of our work that p, = 0. We will compare analytically E {||5^- || 2 } 



and E {||5/i || 2 } in Section VI by evaluating the probability that E |||c)h. || 2 } is greater than E {||<5ft || 2 } when the 



statistics of the channel g are taken into consideration. 

B. Suboptimal Sign Correction 

The sign ambiguity can be resolved suboptimally when only one of the complex channel coefficients, hi, is 
known. The complex channel coefficient he can be expressed in terms of the real coefficients he and hj+e as 
he = he + jhj+e- The suboptimal sign correction is given by 

b s = sgn{h e ue + hj +e u J+ e} = sgn{cos 6 S }, (31) 



where 6 S was used to perform suboptimal phase correction in Section IV We denote by h s = b s u and 5h s = h s — h 
the suboptimally sign-corrected channel estimate, and the resulting error in the estimation of h, respectively. To 
find the resulting MSE, E { ||<5/i s || 2 }, we first write h s in terms of h Q as h s = b s b h . Hence, 

E{||5/i s || 2 } =E{\\b s b K -h\\ 2 } 

= E{||q || 2 }+E{(6A-l) 2 } ( 32 ) 
= E{||q || 2 }+4-4P{6 s = 6 }. 

Using the fact that h e = q ,e + he, and h j+e = Qo,j+e + hj + e, we obtain 

P {b s = b Q } = P {q .ehe + hj + q ,j+ehj +e + h 2 J+e > 0} 
= P {qe + \he\ 2 > 0} , 

where q~e — Qo,ehe + qo,j+ehj+e- The statistics of the random variable q~e are needed to find a closed-form expression 
for P {b s = b a }. It is shown in Appendix [C] that 

q Q ~ j^VrV^Rh, (34) 

\\g\\ 

where V r is a 2J x (2 J — 1) matrix whose columns are the eigenvectors of R corresponding to the eigenvalue 
\, and that 

1 fo 2 ,, „o a 4 



E {q q T } * i^-jUW 2 + T ) V r Vl (35) 
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It can be shown ([34]) that E {q^} ~ and that q # is the sum of N i.i.d. random variables. Using ( [35] ), it is 
straightforward to show that 



1 fa\, ll0 a 4 



E " NW \2 M2 + Tj (N ' " '^ |4) - (36) 
Hence, we may invoke the Central Limit Theorem to approximate q~£ by a Gaussian random variable with the same 
mean and variance, thus obtaining 



I V fllslP + f (i-M' 2 )J 



where Q(-) is the Guassian (^-function [23]. Using ( [37] ), we finally get the MSE expression 



E {«| 2 } - f-rf + ^ (2J - 1) + 4 - 4Q -J 71 JW ) . (38) 

C. Largest-magnitude Sign Correction 



It is clear that the MSE expression in ( |38) decreases as increases. Thus, for a fixed channel g, the MSE is 
lowest when the channel coefficient with the largest magnitude, Hl, is employed in sign correction. Denoting by h a 
and Sh a = h a — h the sign-corrected channel estimate and the resulting estimation error under largest-magnitude 
sign correction, respectively, the resulting MSE is given by 



HI all/ iV|| ff ||4 ^ 2 l|y " 4^ ; ^ Y(f ||flf|| 2 + ^)(l-|^| 2 )/ 

In this case the probability of sign recovery error becomes very small, and the above MSE expression is well 
approximated by the MSE under optimal sign correction in Theorem [2] 

D. Training-based Sign Correction 

In practice, the information needed to resolve the sign ambiguity has to be estimated using training pilots. We 
will derive the probability of making a sign error when we use K pilots to estimate the optimal sign correction b a 
and obtain a closed-form expression for the resulting MSE. We will also provide an expression for the unconditional 
probability of making a sign recovery error taking the channel statistics into account. All K training symbols are 
set to +1. 

The real representation of the K received sample vectors during the training period has the form 

z k = g + n k , (40) 
where k = 1, . . . , K, and g = [R{g} T , Q{g} T ] T . We average the received vectors to obtain 

1 K 

z m = g + — ^2n k . (41) 

k=l 

Using z m , we estimate b Q by b a = sgn{ii T z m } = sgn ^b Q h^ z m |. Hence, 

b = 6 G sgn | h T Q g + ^ ^ hh T n k j 

= 6 G sgn \ \g\ + — h^n k I = 6 D sgn {cos e} , 



(42) 



K 

k=l 



where e is the error in estimating O under training-based phase correction. Denoting by h t = b Q li and 5h t = h t — h 
the channel estimate after sign correction and the corresponding estimation error, respectively, the resulting MSE 
is given by 

¥.{\\5h t \\ 2 } =E{||q || 2 } + 4-4p{6 = 6 }. (43) 



10 



Moreover, we see from d42b that 



p{b o = b o }=pl\\g\\ + ^J2* l onk>0 



(44) 



k=l 



Furthermore, the quantity \\g\\ + -4 Y2k=i ^o n k is a rea l Gaussian random variable with a mean of ||gr|| and a 
variance of which means that 

<2K\\g\W 



P 



Therefore, the MSE in this case is 



E{\\5h t \\ 2 } ~ 



N\\g\ 



2 J 5 r + T )(2J-l) + 4Q 



2^||ff| 



<7" 



(45) 



(46) 



Finally, since the channel coefficients gi,...gj are i.i.d. CN(0, 7 2 ), the unconditional probability of sign recovery 
which takes into account the statistics of the channel g can be found using the Moment Generating Function (MGF) 
of the Gamma distribution [26 ] and is given by 

P u {ho^b } = ^ i 

for K~f 2 /a 2 > 1. 



Kj 2 g (21 



K'j 2 + a 2 



1=0 



(7 



4(K-f 2 + a 2 ] 



(47) 



VI. Performance Comparison 

In this section, we compare theoretically the MSE performance of the two estimators. Our approach is to derive 
closed-forms and/or lower bounds on the unconditional probability that the MSE of the conventional estimator is 
greater than that of the WL estimator, as a function of channel length, SNR and sample size, taking into account 
the statistics of the channel g. We will derive a closed-form expression for this probability for the case of optimal 
phase and sign correction, and two lower bounds on this probability for the case of largest-magnitude phase and 
sign correction. For the other two cases, suboptimal and training-based correction, such probability expressions are 
very difficult to derive. 



A. Relative MSE Performance for Optimal Correction 

Under optimal phase and sign correction, the MSE expressions for the conventional estimator and the WL 
estimator are given in Theorem [T] and Theorem |2| respectively, and the difference in MSE is 

AMSE = E{||J/i || 2 } -E{||«Jh || 2 } 

We are interested in evaluating P{AMSE C > 0}, i.e., the probability that the WL estimator outperforms the 
conventional one. For J > 2, P{AMSE > 0} = P {\\g\\ 2 < a 2 ( J - §)}. The norm-squared of the channel, 
||sr|| 2 , is a Central Chi-squared random variable of order 2 J, which means that J27| 



P{AMSE o >0} = G( ^ (j-l ) ,J\ , (49) 



7 2 V 2 

where G(x,a) is the regularized version of the Lower Incomplete Gamma Function defined in |23| as G(x,s) = 
TTs) fo t s ~ le ~ l dt, where T(-) is the ordinary Gamma function. For fixed s, G(x,s) is strictly increasing with 

respect to x. We will consider the case where ^ > 1, which means that the average received SNR at each antenna 
is greater than or equal to dB. In this case, 

g(^(j- 3 -),j)<g(j-Ij). (50) 



7 2 V 2 J ' J ~ V 2 
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Moreover, it can be shown that G (J — |, J) < ^, which means that the conventional estimator outperforms the 
WL estimator with probability greater than |. At first glance, this is rather surprising because previous works on 
WL subspace-based channel estimation |10j, fl"3| report superior MSE performance for the WL estimator. However, 
none of the previous works considers optimal phase and sign collection. While the WL estimator is indeed superior 
in the remaining three cases (suboptimal, largest-magnitude, and training-based correction), as we will show later in 
this section and in Section |VII| the above result means that the superiority of the WL estimator is highly dependent 
on the imperfect nature of phase and sign correction. The superior performance of the WL estimator under imperfect 
correction makes sense because, intuitively, it is easier to correct the sign which is a binary random variable than 
to correct the phase which is a continuous random variable. On the other hand, when correction of both the phase 
and sign is optimal, the WL estimator is expected to incur more error because the dimension of its observation 
vector is twice that of the conventional estimator. To the best of our knowledge, we are the first to establish the 
relationship between the relative performance of conventional and WL subspace-based estimators and the imperfect 
nature of phase and sign correction. 

B. Relative MSE Performance for Largest-magnitude Correction 

Under largest-magnitude phase and sign correction, the MSE expressions for the conventional estimator and the 
WL estimator are given in ( |22] ) and ([39]), respectively, and the difference in MSE is 

AMSE a = E { ||<5/i a || 2 } - E { ||<5/i a || 2 } 

1 \ a 4 ( J 1 5 N 




4 / (51) 



It is very difficult to obtain a closed-form expression for i-*{AMSE a > 0} using the expression in (51 1. However, 
as we mentioned earlier, the probability of sign recovery error under largest-magnitude sign correction is very small 
and, as a result, E |||<5/i a || 2 } ~ E { ||<5ft || 2 }, which implies that 

AMSE^^fJ^-lU^^ + J^-f) (52) 



N\\g\\ 2 \2\h L \ 2 J iVHsll 4 V 2 2 \ h L\ 2 4 
Using the above approximation, we obtain bounds on P{AMSE a > 0} that are presented in the following theorem. 



Theorem 3 The probability that the WL estimator outperforms the conventional estimator under largest-magnitude 
phase and sign correction can be bounded as follows: 
Case I: (J = 2) 

1 - < P {AMSE a > 0} < 1 - e"^. 

Case II: (J > 3) 

P {AMSE a > 0} > 1 - J 

The proof of Theorem[3]is found in Appendix [Dj We can see from the upper and lower bounds on P { AMSE a > 0} 
for J = 2 that the WL estimator is better at low SNR, while the conventional estimator is better at high SNR. 

For J > 3, we can use the fact that e ~< 2 < 1 to obtain the following lower bound which is looser but does not 
depend on SNR: 

P { AMSE a > 0} > 1 - J ( . (53) 

Both lower bounds that apply for J > 3 approach 1 as J increases for fixed SNR, and the tighter one approaches 
1 as SNR decreases for fixed J. We see from the bound in ((53]) that, starting with J = 4, P {AMSE a > 0} > ~ 
regardless of SNR, which means that the WL estimator performs better with higher probability for this range. For 
J = 3, however, which estimator is better depends on the SNR, as will be confirmed by our simulation results in 
Section |W] 
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VII. Simulation Results 

We use Monte Carlo simulations to verify the accuracy of the analytical expressions we derived in the previous 
sections and to compare the MSE performance of the conventional and WL subspace-based estimators under the 
four scenarios considered in our work. 

We begin with the first three cases: optimal, suboptimal, and largest-magnitude correction. Our MSE results 
are obtained for J = 5 and averaged over the same set of 1000 channel realizations independently generated 
with 7 2 = 1. We show the average MSE performance of the conventional estimator for optimal, suboptimal and 
largest-magnitude phase correction, and that of the WL estimator for optimal and suboptimal sign correction. The 
MSE performance of the WL estimator under largest-magnitude sign correction will not be shown in our plots 
because it is indistinguishable from the performance under optimal sign correction, due to the very low probability 
of making a sign recovery error in this case. In Fig. [2] we plot the average MSE vs. SNR, for N = 100, while in 
Fig. |3j we plot the average MSE vs. N, for an SNR of 10 dB. Figs. [2] and [3] demonstrate that the derived analytical 
MSE expressions are highly accurate. We see that the conventional estimator slightly outperforms the WL estimator 
under optimal phase and sign collection. However, the WL estimator is significantly better than the conventional 
estimator under suboptimal correction and slightly better under largest-magnitude correction. The performance of 
the WL estimator under suboptimal correction approaches the performance under optimal correction for high SNR 
and for large sample size, which does not seem to be the case for the conventional estimator. 

We now consider the practical case where training pilots are used for phase and sign recovery. Using the same 
simulation settings as before, we compare the average MSE performance of the conventional and WL estimator 
when pilot symbols are used to recover the phase and the sign. We study two cases, K = 1 and K = 5. The 
average MSE performance of both estimators vs. SNR under training is shown in Fig. [4j alongside the performance 
under optimal phase and sign correction. As expected, the performance of both estimators improves as the number 
of pilots increase from 1 to 5. However, there is a huge performance gap in favor of the WL estimator. The 
performance of the WL estimator becomes very close to optimal performance as SNR increases for K = 1, and 
becomes indistinguishable from optimal performance for K = 5. The average MSE performance vs. N is shown 
in Fig. [5] for an SNR of 10 dB. In this case, for the WL estimator we only plot the MSE for K = 1 because it 
completely overlaps with the optimal performance. These results show that, in practice, the WL estimator will have 
superior performance because training is much more effective at recovering the sign than recovering the phase. As 
before, both plots demonstrate that the derived analytical expressions are highly accurate. 

In Fig.|6j we plot the analytical expression for P { AMSE D > 0} in ( |49| ) and its experimental evaluation vs. channel 
length for SNR values of 0, 5 and 10 dB. We evaluate P{AMSE Q > 0} experimentally using 10 7 independently 
generated channel realizations (for each value of J) with 7 2 = 1. For each channel realization, we use the analytical 
expressions for the MSE of the two estimators to determine which estimator is superior, since we have already 
established the accuracy of these expressions. We see that under optimal phase and sign correction, the conventional 
estimator is more likely to perform better. As SNR increases, the conventional estimator performs better with 
overwhelmingly high probability. However, as we showed in Fig. [2] and Fig. [3j the difference between the average 
MSE performance of the two estimators is very small. 

Finally, in Fig. |7J we plot the experimental evaluation of P{AMSE a > 0} vs. channel length for SNR values 
of 5, 10 and 15 dB along with the two lower bounds presented in Theorem [5] We evaluate PjAMSEa > 0} 
experimentally using 10 6 independently generated channel realizations with 7 2 = 1. As we did for the case of 
optimal correction, we use the analytical expressions for the MSE to determine which estimator is superior for each 
channel realization. We see from Fig. [7] that, for J = 3, which estimator is more likely to perform better depends 
on SNR, with the WL estimator being favored by low SNR and the conventional estimator being favored by high 
SNR. For J > 4, however, the WL estimator is more likely to perform better regardless of SNR. 

VIII. Conclusions 

In this paper, we presented a thorough theoretical study of the relative MSE performance of conventional and 
WL subspace-based channel estimation in the context of SIMO flat-fading channels employing BPSK modulation. 
Our study explicitly took into consideration the effect of covariance matrix estimation with a finite number of 
received samples and the impact of phase and sign ambiguity resolution. We considered four different scenarios of 
phase and sign ambiguity resolution. The first three assumed that certain information about the actual channel is 
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perfectly known at the receiver. The first one assumed that the phase of the inner-product between the initial channel 
estimate and the actual channel is known, resulting in optimal ambiguity resolution, while the second assumed that 
the phase of a randomly chosen channel coefficient is known and the third one assumed that the phase of the 
channel coefficient with the largest magnitude is known. We derived accurate closed-form MSE expressions for 
each estimator in these cases, in addition to a closed-form expression for the probability that the WL estimator 
outperforms the conventional one in the first case and a lower bound on this probability in the third case. The 
three scenarios resulted in varying relative performances, with the conventional estimator being on average slightly 
better in the first, and the WL estimator being significantly better in the second and slightly better in the last. This 
behavior showed that the relative performance of the two estimators is strongly related to the accuracy of phase and 
sign ambiguity resolution, and that the less information about the actual channel available for ambiguity resolution, 
the larger the performance gap in favor of the WL estimator. We also studied the more realistic scenario where pilot 
symbols are used to resolve the two ambiguities. We derived accurate closed-form expressions for the MSE of the 
two estimators assuming that the same number of pilots are used. Our simulations showed that in this when pilots 
are used there is a significant performance gap in favor of the WL estimator which performs very close to optimal 
even with a single pilot. This scenario showed that the WL estimator significantly outperforms the conventional 
one when the information about the channel which is available for ambiguity resolution is inaccurate. 

Appendix A 
Proof of TheoremQ] 

The phase-corrected channel estimate is h Q = iie 3 ^ u h K Hence, h H h Q = \h H u\ > 0. Since h Q = h+q +a a h, 
we have that h H h Q = l + a , which implies that a is real and a > — 1. Using the fact that h D h Q = 1, we obtain 

a 2 + 2a + \\q f = 0. (54) 

Since a > — 1, the solution of the above quadratic equation is a = + — ||q || 2 - Moreover, since ||q || 2 < 1, 
we can use the Taylor series expansion y/l — \\q \\ 2 = 1 — \ ||<7 || 2 + O (||<Z || 4 ), to obtain a = — §||g || 2 + 
O (Ikoll 4 )' which means that |a G | 2 = j||g || 4 +0 (||<? I| 6 )- Therefore, |a | 2 < ||q || 2 , 5h ~ q a and E {||<5h. || 2 } ~ 
MM 2 }- 

To find E {||q 2 ||}, we denote by A\ the largest eigenvalue of ,R, and let <5Ai = Ai — Ai. Since Rh a = \\h , we 
have 

(R + 6R)(h + q ) = (Ai + 8X x )(h + q Q ). (55) 



Using the fact that Rh = X\h, and ignoring the second-order perturbation terms 5Rq Q and 5X\q in ( [55] ) (as 
in |28|), we obtain 

{R - Xil)q ~ -5Rh + 5\ih. (56) 

Furthermore, we have that (R — X\I) = — ||cjr|| 2 W^, where V is a J x (J - 1) matrix whose columns are the 
orthonormal eigenvectors of R corresponding to the eigenvalue a 2 of multiplicity J — 1 (i.e., the noise subspace). 
Multiplying the two sides of ( [56] ) by W , we obtain 

q Q ~ -^VV H 8Rh. (57) 

He'll 

From eq. (57i, we see thatE{q Q } ~ 0. Let Q = E{q q"}, the MSE is E{||q || 2 } = tr{Q}, and we can expand 
Q as 

Q ~ j^E{VV H 6Rhh H 6RVV H }. (58) 

To obtain a closed-form for the above expectation, we must evaluate the term K{5Rhh H 5R}. Using the approach 
followed in (29j, it can be shown that 

E{5Rhh H 5R} = ^[a 2 \\g\\ 2 hh H + {a 2 \\g\\ 2 + a 4 )/]. (59) 



Therefore, 



Q ^ ^v^p^Hsll 2 + ° A )VV H . (60) 
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From the above equation, we see that, for any £ = 1, . . . , J, E{q ^} = jyjrara (cr 2 ||flr|| 2 + <r 4 )(l — |/ig| 2 ). Finally, 
the MSE is given by 



E {Ikoll 2 } - tr {^p(- 2 ^ll 2 + ° A )W H } = ^(- 2 |^ll 2 + " D- 



(61) 



Appendix B 

In this appendix, we show thatE{g 2 J ~ o for £ = 1, . . . , J. We first let W = E{qq T } ~ E{VV H 5Rhh T 5R*V*Vj}. 
Hence E{q 2 e } = W (£,£), where W(i,j) is the (z, j)th entry of W. To obtain a closed-form expression for W, 
we need to evaluate E{5Rhh T 5R*}. It is easy to see that 

E{8Rhh T 8R*} = ^(E{r{i)r(i) H hh T r(i)*r{i) T } - Rh T R*), (62) 
where r(i) is the ith received vector sample. Moreover, 

E{r(i)r(i) H hh T r(i)*r(i) T } = \\g\\ 4 hh T + 4a 2 \\g fhh T + E{n(i)n(i) H hh T n{i)*n(i) T }. (63) 
To find E{n(i)n(i) H hh T n(i)* n(i) T }, we use the property J30| 

vec(AXB) = (B T ® A)vec(X), (64) 

where the vec(-) operator stacks the column vectors of a matrix below one another and <S> denotes the Kronecker 
product. Hence, we obtain 

E{n(i)n(i) H hh T n(i)*n(i) T } = un\ec[E{n(i)n(i) H ® n(i)n(i) H }vec(hh T )}, (65) 

where unvec(-) converts a vector into matrix format by ordering the elements of the vector into columns of equal 
length. It is straightforward to check that E{n(i)n(i) H ® n(i)n(i) } = c 4 i"j2 + a 4 A where A is a J 2 x J 2 
matrix made of J x J subblocks each of dimension J x J, with the (i,j)th sub-block given by sjsj where Sj is 
the J x 1 vector whose ith entry is 1 and all other entries are zero. Thus, 

E{n(i)n{i) H hh T n(i)*n(i) T } = unvec[(a 4 /j2 + a 4 A)vec(hh T )] = 2a 4 hh T , (66) 

where we have used eq. (39) from (29j. Combining ( [62] ), ( [63] ), and ([66]), we obtain 

l \\g\\ 2 + a 4 )hh T 

(67) 



E{6Rhh T 6R*} = ^(\\g\\ 4 + 4a 2 \\gf + 2a 4 )hh T - ^(\\g\\ 4 + 2a 2 \\g\\ 2 + a 4 )hh q 
= L(2o 2 \\g\\ 2 + <r 4 )hh T . 



Therefore, W ~ -^(2<7 2 ||(/|| 2 + a 4 ) VV H hh T V* V T = 0, since all columns of V are orthogonal to h. Therefore, 

E{<J = W(£,£) ~ 0. ' 

Appendix C 
Proof of Theorem[2] 

Following similar steps to those in Appendix [XJ we can show that \(io\ 2 <C \\q \\ 2 anc ^ that 

Sh ~ q Q ~ ij-i^VrVTiSflh. (68) 

where V r is a 2 J x (2 J — 1) matrix whose columns are the eigenvectors of R coiTesponding to the eigenvalue %r 
of multiplicity 2 J — 1 (i.e., the noise subspace). Thus, E {q } ~ 0. Let 

Q = E {q ql} ~ ^^{V r Vj.E{5Rhh T 5R}V r V^}. (69) 

By applying the same approach used in [29 1 for real signals, we can find the following closed form for E{5Rhh 5R} 
(the derivation has been omitted for brevity): 

E {5Rhh T 5R} = i ((^ 2 ||S|| 2 + hh T + f^-\\gf + l) . (70) 
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Substituting ( [70] ) into ( |69| ), we obtain 

r 2 

iV||g|| 4 V 2"""" 1 ' 4 

Finally, the MSE is given by 



Q-^u lU9\\ 2 + ^r)v r vl (71) 



E {||^ || 2 } ~ tr{Q} ~ ^|| 9 || 2 + £ ) (2J - 1). (72) 

Appendix D 
Proof of Theorem[3] 

The difference in MSE for under largest-magnitude correction is 

AMSE a ^J^('-l) + ^(i + '-^). (73, 



N\\9\\ 2 V2|M 2 J N\\g\\± \2 2\h L \* 4 
Substituting hi = jj~jj"> tne WL estimator outperforms the linear one when 



"fl +-sr^5->^ t-« ■ (74) 



2| 5L |2 11311 2| 5L |2 V4 2 

We first consider the case J = 2, and let p be the index of the channel coefficient with the smallest magnitude 
The condition in (f74]> simplifies to 



2 

M 4 - Isl| 4 + tr 2 |<7 P | 2 + \\9L? > 0. (75) 



It is easy to see from f75| ) that 

P (AMSE a > 0) < 1 - P(\g L \ 2 - \g p \ 2 > a 2 ). (76) 
Moreover, g\ and 52 are i.i.d., so P(L = 1) = P(L = 2) = \ and 

P {\9l\ 2 ~ \ 9p \ 2 > <? 2 ) = \P (M 2 - M 2 > a 2 \L = 1) + ip (| 92 | 2 - M 2 > a 2 |L = 2) 

= P(\ 9 i\ 2 -\ 9 2\ 2 >a 2 \L = l) ( 77 ) 
= 2P(| 5i | 2 -|(72| 2 >ct 2 ). 



Since 51 and gi are i.i.d. CA/"(0,7 2 ), we have that [27| P(|gi| 2 — | ^2 1 2 > <7 2 ) = ?> e r*, which provides us with 
the upper bound 

P(AMSE a > 0) < 1 - e"^. (78) 
Moreover, it can also be inferred from ( f75] > that 

P (AMSE a > 0) > 1 - P (\g\ - \g v \ 2 > ~\ , (79) 

_ ^ 

which allows us to lower-bound P ( AMSE a > 0) by 1 — e ^ . Therefore, 

1 - e"^ < P(AMSE a > 0) < 1 - e~^. (80) 



We will now consider the case J > 3. In this case, the RHS of eq. ( [74] > is negative, and a sufficient condition 
for the WL estimator to be better is 

2\g L \ 2 - \\g\\ 2 < a 2 . (81) 
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It is convenient to define the random variable Z?W = 2\gi\ 2 — \\g\\ 2 , for % = 1, . . . , J. The random variable Z)W is 
the difference of two independent central Chi-square random variables of degrees 2 and 2 J — 2, respectively, and 



P CdM > 



for any x > [27]. Moreover, for any i = 1, . . . , J, 



P(D®>o 3 



P \D® >a 2 \L = i) P(L = i) + P {D® ><t 2 \L^ ij P(L ^ 
P (d® >a 2 \L = i\ P(L = i) 
-p(D®>a 2 \L = 

■\\J-1 



(82) 



Thus, P (£>W > a 2 \L = i) = JP (D® > a 2 ) = J (§) . Going back to the sufficient condition in @T), 

we obtain 

P (AMSE a > 0) > 1 - P (2\g L \ 2 - \\g\\ 2 > a 2 ) 



1 - ^P(2| 9l | 2 - \\g\\ 2 > a 2 \L = i) P (L = i) 



i=l 
J 



l-J2 P ( D ® >a 2 \L = i)P(L 



(83) 



i=i 



J-l 



e i 
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SNR 

Fig. 2. The theoretical and experimental average MSE performance of the two estimators for the cases of optimal, suboptimal, and 
largest-magnitude phase and sign correction plotted vs. SNR for J = 5 and N = 100. 
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Fig. 3. The theoretical and experimental average MSE performance of the two estimators for the cases of optimal, suboptimal, and 
largest-magnitude phase and sign correction plotted vs. N for J = 5 and an SNR of 10 dB. 
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Fig. 4. The theoretical and experimental average MSE performance of the two estimators under training-based phase and sign correction 
using K — 1,5 for the conventional estimator and K = 1 for the WL estimator, plotted together with optimal performance vs. SNR for 
J = 5 and N = 100. 




100 200 300 400 500 

Number of Samples 



Fig. 5. The theoretical and experimental average MSE performance of the two estimators under training-based phase and sign correction 
for K = 1,5, together with optimal performance plotted vs. N for J = 5 and an SNR of 10 dB. 
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Fig. 6. Theoretical and experimental evaluation of the probability that the WL estimator outperforms the conventional one under optimal 
phase and sign correction plotted vs. J for SNR values of 0, 5 and 10 dB. 
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Fig. 7. Experimental evaluation of the probability that the WL estimator outperforms the conventional one under largest-magnitude phase 
and sign correction together with the lower bound in Theorem [3] and the one in eq. < |5 3| > plotted vs. J for SNR values of 5, 10 and 15 dB. 



